arXiv:1504.05584vl [astro-ph.HE] 21 Apr 2015 


Mon. Not. R. Astron. Soc. OOO.IT HTOl f) Printed 23 April 2015 (MN style file v2.2) 

Calcium-Rich Gap Transients: Tidal Detonations of White 
Dwarfs? 

P. H. Sell^*, T. J. Maccarone^, R. Kotak^, C. Knigge^, D. J. Sand^ 

^Department of Physics, Texas Tech University, Lubbock TX 79409, USA 
^School of Physics and Astronomy, Queen’s University Belfast 

^School of Physics and Astronomy, University of Southampton, Highfield, Southampton, SOIL IBJ, UK 


ABSTRACT 

We hypothesize that at least some of the recently discovered class of calcium- 
rich gap transients are tidal detonation events of white dwarfs (WDs) by black holes 
(BHs) or possibly neutron stars. We show that the properties of the calcium-rich 
gap transients agree well with the predictions of the tidal detonation model. Under 
the predictions of this model, we use a follow-up X-ray observation of one of these 
transients, SN 2012hn, to place weak upper limits on the detonator mass of this system 
that include all intermediate-mass BHs (IMBHs). As these transients are preferentially 
in the stellar haloes of galaxies, we discuss the possibility that these transients are tidal 
detonations of WDs caused by random flyby encounters with IMBHs in dwarf galaxies 
or globular clusters. This possibility has been already suggested in the literature but 
without connection to the calcium-rich gap transients. In order for the random flyby 
cross-section to be high enough, these events would have to be occurring inside these 
dense stellar associations. However, there is a lack of evidence for IMBHs in these 
systems, and recent observations have ruled out all but the very faintest dwarf galaxies 
and globular clusters for a few of these transients. Another possibility is that these 
are tidal detonations caused by three-body interactions, where a WD is perturbed 
toward the detonator in isolated multiple star systems. We highlight a number of 
ways this could occur, even in lower-mass systems with stellar-mass BHs or neutron 
stars. Finally, we outline several new observational tests of this scenario, which are 
feasible with current instrumentation. 

Key words: accretion, accretion discs - binaries: close - stars: black holes - white 
dwarfs - galaxies: star clusters - galaxies: dwarf 


1 INTRODUCTION 

The field of time domain astronomy has recently exploded 
with the detection of many new classes of transient sources. 
Considerable progess has so far been made at optical wave¬ 
lengths with the advent of a host of new programs designed 
to search for sources at a variety of cadences. They are also 
more sensitive to finding sources at different luminosities 
in the gap between that of novae and traditional super¬ 
novae llKasIiwalll201ll '). Some of the best-known current tran¬ 
sient surveys are at optical wavelengt hs: the Palomar Tran¬ 
sient Factory fPTF: iRau et al1l2009l l. the Panoramic Sur- 
vey Telescope and R apid Response System (Pan-STARRS; 
iHodapp et al.l l2004l l. and the Public European Southern 
Observatory Spectroscopic Survey for Transient Objects 
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fPESSTO: ISmartt et al.ll2014l '). New observatories (e.g. the 
Large Synoptic Survey Telescope), including some at other 
wavelengths (e.g. the Square Kilometer Array, eRosita), are 
being planned as well. These observatories enable us to bet¬ 
ter map out the parameter space of known transients as well 
as discover new classes of transients. 

One of the recently discovered classes is “calcium-rich 
gap transients” — so called because their spectra contain 
strong forbidden and permitted calcium lines, and they have 
peak absolute magnitudes in the gap between classical novae 
and supernovae. These transients have a set of observational 
characteristics which distinguish them from other classes of 
optical transients. Th ese properties are summarized well by 
llKasliwal et al.ir2012l . hereafter K12): 

• Similar to type I supernovae, showing no hydrogen lines 
in their spectra 
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• Absolute magnitudes in the range from —15.5 to —16.5 
(a factor of about 10 - 30 fainter than the type la super¬ 
novae, which have Mr of about —19.3) 

• Characteristic velocities of approximately 6000 - 
11000 km s“^ 

• Very large calcium abundances, as inferred from their 
nebular spectra 

• Faster evolution than average type la supernovae with 
a rise time ^ 15 days, whereas type la have rise times ~ 
18 days (though it is within the ful l range of type la rise 
times found; e.g. iHavden et al.lliolOh 

• Current, small sample favors the outskirts of known 
galaxies 

As such, many of the properties of these objects are inter¬ 
mediate between those of classical novae and type la super¬ 
novae. 

Previous studies of calcium-rich gap transients have 
used the common relations for type la super novae to esti¬ 
mate some basic parameters of the explosion. IK12| estimate 
the masses of the ejecta by setting the ejecta mass propor¬ 
tional to vt^ (v = photosp h eric v elocity, tr =rise time), as 
first suggested by lArne 3 ||1982D , and scaling to type la 
supernovae. This yields typical total masses of ~ 0.3 - 
0.7 Mq, where the fraction of ®®Ni produce d is ver y small 
(e.g. 0.016 Mq for PTF lOiuv). However, |A meld ’s mod¬ 
elling explicity requires that there be substantial ®®Ni in 
the supernova ejecta, and modelling of the nebular spec¬ 
trum o f SN 20 05E in |K12| gives a mass about half as large 
as the lArnetd formula does. Type la supernovae are pre¬ 
dominantly powered by radioactive decay of iron-peak el¬ 
ements, where the amount of nickel produced would be 
large and the amount of c alcium produced wo uld be small 
(e.g. IWooslev et al.l [l98^ : iBildsten et al.l[2007h . Therefore, 
the modelling falls short in properly describing the explo¬ 
sion. 

Various models exist to explain the unusual character¬ 
ises of these transients, many of which are summarized in 
Hi Many of them have at least one serious flaw that makes 
them unable to explain these transie nts. To explain the 
calcium-rich gap transient, SN 2005E, IPerets et al.l (120101 1 
concluded that they were witnessing helium detonation in an 
interacting double white dwarf (WD ) system with a helium 
WD mass donor llShen fc BildsterJ l2009l : IWaldman et al.l 
l20ldl . but the light curve is not well-matched. Various chan¬ 
nels fomthndetonation_ofs^-Chnndrasekhnrmmss_V^s ex¬ 
ist (IWooslev fc Wea^^ll994 IWooslev fc KasenI 1201 ill , but 
such models do not provide a good match to the observed 
light curves. Accretion-induced collapse (AI Cl of a rapidly 
rotat ing white dwarf into a neutron star (iMetzger et al.l 
I 2 OO 9 II is also ruled out based on the predicted shape of the 
light curves, as well as the large velocities an d abundances o f 
intermediate elements predicted. Recently, iMetzgeJ (l2012l l 
presented a one dimensional nuclear-dominated thick accre¬ 
tion flow (nuDAF) model resulting from the unstable mass 
transfer following Roche lobe overflow, but it does not pro¬ 
duce the right amount of helium and calcium to match many 
of these transients. Some of the difficulty in finding a “best- 
fit” model may arise in the heterogeneous nature of the cur¬ 
rent sample: there may be selection bias in finding these 
transients far from their host galaxies, a range of velocity 


widths in the lines in the nebular spectra, different light 
curve e volutions , for e xamples. 

In iMetzgerl (l2012ll . another scenario is mentioned: the 
tidal disruption of a WD by a neutron star or black hole 
(BH) leading to a nuclear runaway (detonation). However, 
this scenario is not discussed further. Some other work 
has considered the tidal disruption of a WD in other con- 
texts but without connectio n to calcium-rich gap transients. 
IClausen fc EracleousI (l201lll considered the emission spec¬ 
trum after the disruption but not detonation of a WD by a 
BH; the spectrum is incomplete because they did not con¬ 
sider the fusion of heavy elements in the tidal debris (con¬ 
sidered only mass fractions of 67% oxygen, 32% carbon, and 
1% helium). While inves tigating the vertica l compression of 
tidally disrupted debris, IStone et al.l (l2013lj noted the pos¬ 
sible importance of WD tidal disruptions to gr avitational 
wave signals. Finally, IShcherbakov et ^ (l2013l l proposed 
that a highly unusual combination of a gamma ray burst 
and a supernova explosion could be explained by a WD tidal 
disruption. 

In this work we consider the tidal disruption leading to 
nuclear runaway (detonation) of WDs in more detail to ex¬ 
plain the calcium-rich gap transients. The goal of the paper 
is mainly to discuss the plausibility of this scenario in connec¬ 
tion to calcium-rich gap transients, something that has not 
yet been done in the literature. Future more detailed analy¬ 
ses of this scenario in the context of Ca-rich gap transients 
will need to be undertaken when better multiwavelength (es¬ 
pecially high-energy) observations become available. These 
are needed to solve some of the large uncertainties and de¬ 
generacies in various model parameters discussed through¬ 
out the text. 

In section [2] we describe the model and how it can 
explain various observational characteristics: masses, abun¬ 
dances, light curves, and the positions of these sources. Then 
in section [3] we discuss some possible scenarios that could 
produce the tidal detonation and calculate expected rates of 
these transients to verify if the proposed origin is plausible. 
In section [4] we then analyze a Chandra observation of one 
of these transients, which enables us to place upper limits 
on the masses of the WD and BH in the progenitor system. 
Finally in section (5] we summarize our results and suggest 
observational tests of our model. 


2 MATCHING THEORY WITH 
OBSERVATION 

We p ropose that the model developed by iRosswog et al.l 
ll2009l . hereafter R09) naturally explains the set of calcium- 
rich gap transients. R09 use a three-dimensional smoothed 
particle hydrodynamic (SPH) method to simulate the close 
interaction between a WD and a more compact object (a 
BH) that completely disrupts the WD. When a WD passes 
too close to a BH (within the Roche or tidal disruption ra¬ 
dius), the differential gravitational force from the BH will 
overcome the gravitational force of the WD. This tidal force 
can elongate the WD in the orbital direction and compress 
the WD in the direction perpendicular to the orbit to a high- 
enough density and pressure to initiate runaway nuclear 
burning. If the nuclear energy is greater than the binding 
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energy of the star, an explosive nuclear detonation occurs, 
producing a luminous transient eveirt. 

The tidal disruption leading to the possible detonation 
of a WD will only occur when Mbh < 2 x 10® Mq or 
stellar- or intermediate-mass Bh£], (IMBHs; see R09 for a 
more detailed discussion). R09 consider 100 Mq < Mbh < 
10000 Mq, but BHs at a larger range of masses (up to 
Mbh ~ 10® Mq or Mbh < 100 Mq) or even neutron stars 
may be viable detonators as well. Next, we discuss how each 
of a set of testable theoretical predictions of this model can 
be explained by current observations of this group of tran¬ 
sients. 


2.1 The mass of the detonated white dwarf and 
abundances of leftover material 

There already exist a set of strong theoretical and observa¬ 
tional ways to constrain the mass of the detonated WD in 
this scenario. We first explore the upper limit to the mass of 
the WD, which is most strongly motivated by the abundance 
patterns in the nebular spectra. 

As implied by the name, these calcium-rich gap tran¬ 
sients produce a relatively large amount of calcium. How¬ 
ever, they do not produce a lar ge amount of iron-peak ele¬ 
ments (see the spectra in lK12h . This is only well-m atche d 
to the t idal detonation scenario if Mwd < 0.6 Mq ilRogh . 
As IROQI present the most realistic description of a tidal det¬ 
onation of a WD to date, this is currently the best estimate 
of reaction rates in this scenario. However, there are some 
important limitations to this work. The upper limit to the 
WD mass is not very precise because the WD mass grids 
presented in are very coarse and might be a function of 
other parameters (e.g. tidal radius relative to pericenter, BH 
mass). In addition, the exact boundary of the upper limit is 
also not well-determ ined because the nucleosynthesis calcu¬ 
lations presented in IrooI are based on a simplified nuclear 
reaction network. 

This latter point is particularly important because this 
implies one cannot make specific abundance predictions 
from this model. For instance, silicon is taken to represent all 
the elements within the quasi-equilibrium group near silicon 
in this network. Therefore, given that calcium lies within 
this group, as does scandium, it should not be a surprise 
that large amounts of calcium are not specifically predicted. 
The tenta tive detection of st rontium in one calcium-rich 
transient (ISullivan et ahlBoilIl would potentially present a 
problem, as strontium cannot be produced in this scenario, 
but it should be noted that this detection is of marginal 
significance and only has been seen in one object. A clear 
prediction of this picture is that, as a sample of events with 
well-measured abundances emerges, progressively higher ra¬ 
tios of the iron-peak elements to the elements near silicon 
in atomic number will be seen as the peak luminosities rise. 
However, since we expect a range of impact parameters for 
the tidal detonation model, future work may require a very 
large sample size to discern these trends. 

In addition to the abundance analysis, which requires 


^ WDs are not tidally disrupted by supermassive BHs because 
the tidal disruption radius is inside the event horizon; they are 
swallowed whole. 


that Mwd < 0.6 Mq, the density-dependent Roche limit 
may favor the tidal disruption of lower-mass WDs. Lower- 
mass WDs have lower densities (p) because they are larger. 
Because the minimum distance at which an object is dis¬ 
rupted by tidal forces scales as this may imply that 

lower-mass WDs are more susceptible to stretching and com¬ 
pression, possibly making them easier to deto nate. However, 
based on the limited number of grids in IROQI . it is not clear 
whether this susceptibility outweighs other effects (e.g. need¬ 
ing to reach a critical density to detonate, which would re¬ 
quire more compression); more work on this topic is clearly 
needed. 

If this is true, then this would imply that the lowest- 
mass WDs are be preferred in this scenario. This is a prob¬ 
lem for single star evolution, as it only produces WDs 
with M > 0.45 Mq for the lowest-mass evolved stars i n 
the universe (e.g. iMarsh et al.l 1 19951 : I Kepler et al.l l2007l l. 
However, binary processes can produce lower-mass systems 
(common envelope ev olution and Roche lobe overflow; e.g. 
iNelemans et al.ll200ll ~) and account f or the low-mass WDs 
found in the last cou ple decades (e.g. iBergeron et al.l [19921 : 
iBragaglia et aI.lll995l L The lowest-mass WDs yet found are 
commonly referred to as extremely low mass (ELM) WDs 
~ 0. 2 Mq and are he l ium-dominated (e.g . jLiebert_^^_^ 
2004; iKilic et al.l l2007l : IVennes et al.l l201lf : iHermes et al. 
2013l) . Putting these constraints together, we conclude that 
0.2 Mq ^ Mwd ^ 0.6Mq. 


2.2 The continuing power source and light curve 


Absorption lines in type la SNe tend to become narrower 
with time, as the outer layers become more diluted by the 
expansion of the ejecta and the line-forming region recede s 
to lower velocities (e.g. lFilippen^ll997l : iFolev etld] 1201 il l . 
However. I Valenti et al.l (2014) observe the opposite trend in 
one of the calcium-rich gap transients, SN 2012hn: constant 
photospheric velocity (as inferred from the blueshift of ab¬ 
sorption lines) and increasing linewidths, which are better 
explained by a continuing energy source. 

In supernovae, the radiative decay of nickel serves as 
a continuing source of energy after the original ejection of 
material in the explosion. However, as discussed in the previ¬ 
ous section, iron-peak elements are not appreciably seen for 
calcium- rich g ap transients nor would we expect them in this 
scenario. IR09| find that, in the low-mass WDs (~ 0.2 Mq), 
only about 0.03 Mq of iron-peak elements are likely to be 
produced in the explosion, and the detonated 0.6 Mq model 
produces an abundance of iron-peak elements ~ 4x less than 
the detonated 1.2 Mq models. 

Instead, the power source can be naturally explained 
by fail-back accretion of a fraction of the WD material onto 
the BH. The very hot inflow of gas (a few tens to hundreds 
of million K) not only produces an optical flare but also 
produces copius X-ray emission. The highly blueshifted and 
broadened lines (up to thousands of km s“^) can be ex¬ 
plained by motions of material associated with the accretion 
flow, as seen for both s tellar and supermassive BHs (e.g. disk 
wind s; Murray et alJ 19951 : iMiller et al.l l2008l : iKing et ^ 
l2012l:lTombesi et al.ll2015lL 

For instance, IroqT show that for the detonation of a 
0.2 Mq WD by a 1000 Mq BH, the accretion rate onto the 
central BH will initially be very large (~ 10~® "^ Mq s“^) 
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for the first approximately 10 min and then decay as 
dReedllO^ : IPhinnevUlQSgl ). Detailed modelling of the opti¬ 
cal light curve using type la templates has been challeng¬ 
ing for all previous work and is beyond the scope of this 
paper. Nevertheless, a comparison of the slope of late-time 
R-band light curve 10 - 100 days after the peak) of a few 
of the calcium-rich gap transients with well-sampled light 
curves (e.g. SN 2005E, RTF lOiuv) indicates that the tidal 
detonation scenario (M oc seems to provide a bet¬ 

ter matc h tha n the exponential radioactive decay of ®®Ni as 
done in |kT3, even when taking into account more of the 
decay chai n (^^Co and the ki netic energy loss from its (3+ 
decay; e.g. IPado fc Darj|2oi^ '). Another more more compli¬ 
cated possibility is that the optical light curves could be 
explained by the decay of lower-mass elements (^®Cr — 

—^'*®Ti and ^^Ti —^-^^Sc —^^'*Ca), producing late-time emis¬ 
sion dValenti et al.llioi^ l. 

Later on, when the accretion rate falls be¬ 
low the Eddington rate onto this BH, Lsdd ~ 
10^^(MB/f/10® M©) erg s“^ or Msdd = 2.2 x 

10 “®(M_bb/10® Mq) M f:] yr~^ (standard radiative ef¬ 
ficiency 0.1 assumed; e.g. Novikov fc Thorng|l973l L which 
occurs approximately seven months after the deton ation 
of a 0.2 Mq WD around a 1000 Mq BH in IROOI . the 
X-ray luminosity should begin to decay. However, the exact 
brightness and evolution of the X-ray bright phase depends 
a number of poorly constrained parameters: the WD’s 
initial mass and trajectory, the mass of the BH, and the 
stru cture of the initially super-Eddington accretion flow 
fe.g. iLodato fc Rossill201ll : ljiang et al.ll2014l ). for examples. 

As an illustration of the importance of the assumed 
structure of the accretion flow, the viscous accretion 
timescale is proportional to the rati o of the disk sc a le height 
to tid al disruption radius squared dLi et al.| [ ioo^; iMetzgerj 
I2OI2II . Accretion disks can radically vary in scale height de¬ 
pending on the accretion mode, which is tied to the effi¬ 
ciency of radiative cooling (iHavasaki et al.l[2015h . These un¬ 
certainties make the viscous dissipation timescale uncertain 
by many orders of magnitude. This strongly affects our abil¬ 
ity to apply upper limits to physically meaningful parame¬ 
ters {Mwd, Mbh) in our X-ray observations of one of the 
calcium-rich gap transients in section [4] 


Mr > —5. 6 mag for SN 2012hn, and Mr > —7.3 mag for 
SN 2003H; I Lyman et al .112013 1^ 

For these systems, these observations make the sce¬ 
nario of the tidal detonation of a WD caused by a random, 
flyby encounter with an IMBH very unlikely because the en¬ 
counter rate would be too small without the presence of a 
dense stellar system. In these cases, interactions pro ducing 
high- velocity runaway systems must be considered jFolevI 
l2015ll . where the detonations occur in isolated three-body 
systems. We provide a detailed discussion of random en¬ 
counter rates and possible three-body interaction scenarios 
that can satisfy these requirements in section [31 


2.3.1 Can these stellar systems survive? 

To discern whether globular clusters and dwarf galaxies will 
survive long enough to be observable in the local universe, 
we can apply the Roche or tidal disruption limit to these pos¬ 
sible host stellar systems. The dynamical friction tim escale 
for galactic cannibalism is jBinnev fc Tremainell200^ ') 


tfric — 


2.7 Gyr n 


/ OM W 1 

^100 km s"^\ 

V200 kms-i/ ' 

V CTs / 


( 1 ) 


where A is the Coulomb logarithm, ri is the initial radius, 
(Tm is the typical velocity of these hosts through the halo, 
and (Js is the veloci ty dispersion with i n a gl obular cluster or 
dwarf galaxy. If the lFaber fc Jacksorj (Il976l~) law is satished, 
then the value of A will be 2^^'^gm/ us I Binnev fc Tremaine! 
l2008lfl . We pick an intermediate velocity dispersion between 
the largest dwarfs and globular clusters, ag — 40 km s“^. 
Using the typical di stanc e of these objects from their host 
galaxy (n ~ 25 kpc. lKljf , we find that it would take about 
a Hubble time for these hosts to inspiral if they start on 
a circular orbit (and eccentric orbits would only have semi¬ 
major axes larger than this for the objects to have dynamical 
friction timescales shorter than a Hubble time). Therefore, 
many of these objects will exist in the semi-local universe. 
Because globular clusters typically have higher central stel¬ 
lar concentrations, they will surviv e to much smaller g alacto- 
centric radii than dwarf galaxies (iMaver et aI.|[200a L how¬ 
ever, uftrar^comg^tdwarf galaxies might be an exception 
(e.g. IStrader et al.ll20i^ . 


2.3 Positional information 

To date, the calcium-rich gap transients have been located 
predominantly on the outskirts of their host galaxies (if this 
is not a sele c tion b ias, which is not clear at the present time). 
[Yuan et al.l (l2013l 'l find that the transients are usually well- 
beyond the extent of the nearby host galaxies (-^ 5x the 
K-band half-light radius on average). They conclude that 
these transients do not trace the stellar mass profile but 
could be matched to globular clusters or dwarf galaxies. 

Initial non-detection limits for many of the sources 
{Mr > —10 - 12 mag) required that h ost stellar clusters 
have L < 10^ Lq (iKasIiwal et al.|[2012lL However, recent, 
much deeper observations searching for hosts for a couple of 
the transients have ruled out all but the faintest dwarf galax¬ 
ies and globular clusters {Mr > —5.3 mag for SN 2005E, 


3 SPECIFIC ORBITAL SCENARIOS AND 
RATES DISCUSSED 

As discussed in the previous section, likely host systems for 
these transients are faint dwarf galaxies or globular clusters 
given our scenario. Precise rate predictions for this scenario 
cannot be calculated at the present time due to consider¬ 
able uncertainties in multiple key parameters: the fraction 

^ SN 2003H’s upper limit is not as deep because it lies in a region 
of higher background between two merging galaxies. 

® The fainter end of the dwarf galaxy luminosity function shows 
some evidence for increased dark matter fractions and mild de- 
viations from the fu ndamental plane for elliptical galaxies (e.g. 
IZaritskv et al.ll2006h . but this should not affect our qualitative 
conclusions. 
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of globular clusters and faint dwarf galaxies which contain 
intermediate- or stellar-mass BHs, the masses of those BHs, 
the masses of the WDs that might be detonated, the struc¬ 
tural parameters (and hence central stellar densities) of the 
clusters and/or faint dwarf galaxies which contain those 
BHs, etc. We only consider rough estimates to test if this 
model is plausible or not. A thorough analysis is well be¬ 
yond the scope of this paper and is very challenging because 
it would have to be quite complex, taking into account (not 
an exhaustive list): the binary fraction around the primary 
BH, general relativistic corrections, stellar densities, the re¬ 
plenishment rate of the loss cone of the BH, the spin of the 
BH, the mass ratio of the BH in the case of a binary, etc., 
while having good observational constraints for the relevant 
parameters. 

We consider two basic plausible scenarios in which a 
WD could get close enough to a BH with Mbh iS 10® Mq 
to be tidally disrupted: 1) the chance flyby of a single WD 
to an IMBH and 2) the perturbation of a three-body sys¬ 
tem, which is especially needed for the stellar-mass BHs with 
lower interaction cross-sections (oc R%h oc M%b)- Then 
we consider how these ideas fit with the locations of these 
events. 


3.1 Chance white dwarf flyby of an 
intermediate-mass black hole 

Given the larger interaction cross section of an IMBH, 
chance flybys in dense stellar systems could occur often 
enough to be observable calcium-rich gap transients. Rate 
estimates testing for plausibility are easier to make for glob¬ 
ular clusters, so we consider them first, despite the fact that 
the positions of the transients far from the centers of bright 
nearby galaxies may favor faint dwarf galaxies as the sites 
of the explosions. 


3.1.1 Rates in globular clusters 

iBaumgardt et aP (120041 ') estimates a rate of ~ 10“® WD 
tidal disruptions per globular cluster per year, assuming 
typical clusters have a 1000 M© BH and a central stellar 
density of 10® stars parsec“®. Given ~ 400 globular clus¬ 
ters per typical Milky Way-like galaxy (taking a number 
larger than that for the Milky Way, since it s specific fre¬ 
quen cy of globular clusters is relatively low, iGalleti et al.l 
[2003), and a density of galaxies of about 1 Mpc“®, then 
within 200 Mpc, there should be ~ 10^° globular clusters 
within the distance range out to which the calcium-rich gap 
transients have been seen. If ~ 10% contain IMBHs, then 
it would be reasonable to expect ^ 1 tidal detonation per 
year. Given that t he tidal disruption rate goes as the square 
of the BH mass llBaumgardt et all l2004l l , the rate would 
be expected to increase in the event that there is a broad 
distribution of BH masses. The crude rate calculation can 
be taken as evidence of plausibility, under the assumption 
that some non-negligible fraction of globular clusters contain 
IMBHs. 

At the present time, the evidence that globular clusters 
contain IMBHs is weak. Some claims have been made on 
the basis of stellar dynamical evidence for an increase in the 
charateristic velocities of stars in the centers of a few clusters 


fe.e.lNewell et al.lll976l: 

Gerssen et al.ll2002l: iGebhardt et al.l 

I 2 OO 2 I: Novola et al. 2008 

). On the other hand, in most cases. 


the data can be equally well explained by mass segregation 
of neutron stars and WDs down to the central regions 
of the cluster, causing an incr ease in the mass-to-Iight 
ratio without an IMBH (e.g. Illlingworth fc Kind Il977l : 
IBaumgardt et aP l2003l 'l. Alternatives to radial velocity 
searche s, such as proper motion searches for dynamical evi¬ 
dence llMcLaughlin et ^l2006l : lAnderson &^_TOi 2 _der MA^ 
2 OI 0 I) and searches fo r radio emission dMaccaxone 


2004 


Maerau;one_^^_^ 2005^ de_^i]cke_etalj 2006: Bash et all 


20081: Maccarone fc Servillat I 2 OO 8 I: ICseh et al.l I 2 OI 0 I : 

Lu fc Kongl 20lil ~ Strader et al.l 2012ll have failed to find 


any IMBHs in Galactic globular clusters. The M31 cluster 
G1 does s how both dynamical evidence for a BH and radio 
emission (lUlvestad et al.|[200itl . although more recent radio 
follow-ups have failed to detect a source at that position, 
indicating either that the detection was spurious or the 
source is strongly variable dMiller-Jones et aDl2012h . 

Despite the lack of evidence thus far for IMBHs 
in globular clusters, IMBHs may have been more easily 
formed through the formation of more massive stars i n 
lower metallicity environments (^e.g. iBelczvnski et akllioi^ l. 
This may also be helpful because it does not require 
the merger of stellar-mass BHs (e.g. iMiller fc HamiltonI 
l2002ll. where the grav itational radiation rocket effect 
( RndmounG^Jleej[l989|) can eject t he IMBH from the clus- 
ter (iHolIev-Bockelmann et al.l[200^ 1. Indeed, deep observa¬ 
tions of globular clusters around the Galaxy and other galax¬ 
ies over the last few decades have discovered globular clus¬ 
ters at large galactoce ntric radii and shown that they have 
low metallicities (e.g. ISearle fc ZimJ Il978l : iRhode fc Zeoll 
I2n04h . However, whet her or not t he globular cluster radial 
composition gradient dZin nlll985lf ex tends farther out into 
the halo is not clear dVan Den Bergh|[2003tl . 

As discussed Section [2Hl the tidally disrupted object 
must be a low-mass WD and is likely a helium WD. Un¬ 
fortunately, the number density of helium WDs in globular 
clusters is not well-known. Helium WDs cannot form at the 
present epoch through normal single-star evolution; they can 
only form through binary evolutionary processes (see e.g. 
iHansen et al.l 1^031 ). In globular clusters, this can happen 
through a direct collision between a red giant and a compact 
star. This mechanism has been proposed for producing the 


ultra compact X-ray binaries in globular cli:isters_JVerbi^ 


19871) and is probably at work in NGC 6397 (|Grir]dIa^_et_^ 
200 ll : I Taylor et al.ll20(5H : iHansen et al.ll2003l : IStrickler et al 
2OO9II . 

Placing observational constraints on the density of he¬ 
lium WDs in nearby globular clusters is extremely difficult 
because they are so faint and may not be possible because 
they should escape from the centers due to mass segregation 
(e.g. iFregeau et al.l[200^ . Only in the presence of an IMBH, 
which appears to be absent in these systems, is mass seg¬ 
regation quenched because the interactions with the IMBH 
domi nate over the interactions with other stars dGill et al.l 

I 2 OO 8 II . 


3.1.2 Rates in dwarf galaxies 

Observations clearly show a dearth of dwarf galaxies within 
a few tens of kpc and those found in this region are 
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highly disturbed (e.g. Saggitarius Dwarf, M32, Willman 1 


Ursa Major I I; BekMeLah 200l|; Johnston et al.l 1 19951 : 


IWillman et alJ l2005l : iMunoz et ^ l2010l i . The “classical” 
dwarf spheroidal galaxies orbiting the Milky Way are all 
located at least 25 kp c from the Galactic center (e.g. 
iGallagher fc Wv^ll994l. as are t he new ultra-faint dwarf 
galaxies (e.g. Belokurov et al M) . Streams of past disrup¬ 
tion events of dwarf galaxies by the Milky Way clearly sug - 
gest this is still occurring locally (e.g. Ijohnston et al.lll99a) . 
The preferential locations of faint dwarf galaxies are in con¬ 
trast to the distribution of globular clusters, which includes 
many objects within a few kpc of the Galactic center. There¬ 
fore, the finding that the few calcium-rich gap transients so 
far are located at large galactocentric radii provides early 
suggestive evidence in favor of faint dwarf galaxies rather 
than globular clusters as the host sites (barring a selection 
bias). 

Thus far, the best case for an IMBH in a small dwarf 
galax y is that of HLX-1 near ESO 243-49 (iFarrell et al.l 
l2009l j. This system shows X-ray emission that peaks about 
10^^ ergs/sec , show s state transitions like stellar mass BHs 
dGodet et al.llioogll . is spectroscopically confirmed to be a 
membe r of the cluster of ga laxies against which it is pro¬ 
jected (IWiersema et al.ll201(]l l. and recently has been found 
to be in a cluster or faint dwarf gal axy which likely has a 
young population (IFarrell et al.ll201^ j. It has been suggested 
that the companion star to the possible IMBH may have 
been tidally captur ed, given its rather long orbital period 
llLasota et al.lBoill I. 

Next, we proceed to test this scenario for plausibility as 
we did for globular clusters, although the rate estimate for 
tidal detonations is quite crude because the input parame¬ 
ters are so poorly constrained. We use the tidal disruption 
rate relation from iBaumgardt et akl ll2004h , where the key, 
variable parameters are the IMBH mass (Mimbh), the cen¬ 
tral stellar density (ric), and the central velocity dispersion 
(Uc). Given an upper limit on the m ass of the host galaxies of 
many of these even ts of 10^ Mp) and the BH-bulge 

mass relation (e.g. iMcConnell &: mJ 2013ll . dwarf galaxies 
would have Mimbh ~ a few x 10^ M©, roughly consistent 
with the upper limit derived in Section (4] A tidal detona¬ 
tion would occur in dwar f galaxies with dense central stellar 
clusters like M60-UCD1 (IStrader et al.l[2013l l but at least a 
factor of a few fainter. Considering that the central surface 
density of M60-UCD1 is close to but slightly less than that of 
the densest globular clusters, we adopt Uc ~ 10^ stars pc“®. 
We adopt Vc = 40 km s~^ as in section assuming that 
a is not considerably different from Vc- 

Assuming the t idal disruption rate ro ughly scales as 
(IBaumgardt et all 120041 '). we find that 
the dwarf galaxies with IMBHs have a tidal disruption 
rate per object about 20 times higher than globular clus¬ 
ters. Then, assuming 1/10 globular clusters have IMBHs, 
we would need only 1/200 as many dense dwarf galaxies as 
globular clusters to produce the observed rate of calcium- 
rich gap tran sients. Surveys of the se galaxies are far from 
complete but IPhillipps et akl (120011 1. for example, estimate 
that there are ~ 100 dwarf galaxies in the central region of 
the Fornax cluster in the magnitude range where they would 
be plausible hosts for IMBHs, and this region contains con¬ 
sider ably less than 20000 globular clusters (e.g. lOstrov et all 
Il998l l. Thus, faint, dense dwarf galaxies appear to be a vi¬ 


able candidate site for tidal detonations of helium WDs by 
IMBHs for at least some of the calcium-rich gap transients. 

3.2 Three-body interactions 

There are a variety of ways in which three body interactions 
could perturb a WD to send it into a close encounter with 
a B H. Such interactions are frequently intrinsically chaotic 
(e.g. IPortegies Zwart fc Boekholf 20141'). even resulting in 
orbital flips, for instance ( Li et alT 2014^ . Given the very 
larg e body of literature regarding three-body interactions 
(see lValtonen fc Karttunenll2006l : iMusielak fc Quarlesll2014l . 
for summaries), we only briefly describe and discuss a few 
possible scenarios that would induce the tidal disruption of 
a WD. Again, this summary is meant only as a demonstra¬ 
tion of plausibility for three-body interactions as a possible 
channel for the WD tidal detonations, not an exhaustive 
discussion or a detailed analysis of any particular model. 

While the study of the three body problem goes 
back hundreds to years to as early as the days of Sir 
Isaac Newton, considerable development on the problem 
has occurred in approxim ately the last half century (e.g. 
lEggleton fc Kiseleval Il995lf . Much work has been focused 
on the well-known Kozai-Lidov mechanism, where perturba¬ 
tions in the three-body system transfer angular momentum 
betwee n orbital modes with dif ferent eccentricity and incli¬ 
nation (iKoza ill 19621: lLidovlll962lI . Since then, numerous sim¬ 
ulations (e.g. Anosovall 19901 1 have been carried out to explore 
the large fa mily of solutions tha t exist, from larger N-body 
encounter s (iLeigh fc Geli^l2012l l to higher order expansions 
(octupole. iNaoz et al.l 2013l l. The importance of three-body 
interactions has also been explored for a variety of astro- 
physical contexts outsid e of the solar syst em (e.g. ultra- 
luminous X-ray sources, iBlecha e t al.l l20 06l: blue str aggler 
stars in dense gl obular clusters, CTiutteriee et al ]|2013l: grav- 
itational waves, lAntonini et aklboiJl : the final parsec prob¬ 
lem involving the merger of supermassive BHs, iKhan et akl 


2 OI 3 I: the motions of stars ar ound BHs in galax y centers, 
Sesana et al.ll20()^ . l2007l . 1200^: extrasolar p lanets, iMardlina 


2 OI 3 I : white dwarf mergers, Hamers et ahll^lSl l. We now 


consider its importance for the tidal disruption of WDs when 
they are paired with a BH and another member (another 
star or remnant) in a triple system. We consider paths to 
exchange angular momentum which do and do not involve 
mass transfer. 

First, we consider a few examples of three-body in¬ 
teractions that will exchange angular momentum without 
mass transfer. The first two options require the Kozai- 
Lidov mechanism and are enhanced i n higher-mass-ratio 
systems. Rapid eccentricity oscillations Jlvanov et al.llioosl 'l 
can make the magnitude of angular momentum in the in¬ 
ner binary comparable to the total angular momentum of 
the system. Wh e n num erically evolving the triple system, 
lAntognini et al.l (I20l4l find that, because the general rel¬ 
ativistic terms are sensitive to the orbital eccentricity, the 
inner binary can merge well within a few Gyr, or even faster 
if the mass ratio is large because it is proportional to the am¬ 
plitude of the eccentricity oscillations; this would be the case 
for a WD and BH inner binary. Second, a BH binary system 
in a stellar cusp strongly enhances the rate of t idal disrup¬ 
tion over t wo-body relaxation for a single BH dChen et akl 
l2009l . I^llll . This was originally demonstrated for a super- 
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massive BH and another BH in a binary in the center of a 
stellar cusp of a galaxy, but also applies to a stellar-mass BH 
and IMBH pair in the center of a globular cluster or dwarf 
galaxy, which are relevant for calcium-rich gap transients. 
The tidal disruption rate in this case is at least approxi¬ 
mately an order of magnitu de larger than a single IMBH 
fed by two-body relaxation dWang Sz Merritt l2004h but it 
also quickly depletes the loss cone in < 10® yr. 

Yet another way a triple system can cause a WD in 
the inner binary to migrate toward a BH is through the 
tidal excitation of high frequency p-modes. The se can redis¬ 
tribu te angular momentum in a triple system dFuller et al.l 
I2OI3II . causing the inner binary to inspiral. Finally, the re¬ 
distribution of angular momentum does not have to occur 
within a triple star system to significantly perturb a close 
binary. A close flyby of another star/remnant or binary can 
induce a high-enough eccentricity in a WD-BH binary that 
forces the WD to pass within the tidal d isruption radius (see 
iHeggie fc Rasid fToOd : Im accaron j l2005l . for detailed discus¬ 
sions). 

Finally, we consider two exa mples which in volve mass 
transfer. In the first of these cases. I^ahleil (120101 1 show that 
the mass transfer/accretion of gas from a dense interstellar 
medium (star-forming region) can act as a way to transfer 
angular momentum. When moving through and accreting 
the dense gas, the orbit can be shrunk by dissipating angu¬ 
lar momentum through acoustic waves. This would require 
the site of the transient to be near a dense molecular cloud; 
however, these transients are not g enerally associated with 
star formation (iLvman et al.ll^lSll . Another possibility is 
that mass is transferred from star to star within the bi- 
nary or lost as the star s evolve off of the main sequence 
dMichaelv fc Peret j|2014l l. This can happen in a much older 
stellar population and induce a chaotic orbital behavior. 

As has just been made clear, there are a multitude of 
various ways three-body interactions can lead to the tidal 
disruption of a WD, many of which do not require extremely 
dense stellar systems (e.g. the importance of three-body 
inter actions in lower-density stellar groups; iLeigh fc Gellerj 
I2OI3II . This makes it very difficult if not impossible to rule 
out the tidal disruption origin on interaction cross-section 
arguments alone. 


4 THE PARTICULAR CASE OE SN 2012HN 

To test the proposed origin of the calcium-rich gap tran¬ 
sients as tidal disrup tions of WDs, we ob served one of these 
sources, SN 2012hn d Valenti et al.l[2014l l to look for accre¬ 
tion emission during the decay phase. SN 2012hn has all of 
the typical characteristics of calcium-rich gap transients as 
outlined in Section [T] It is also fairly nearby, at a distance 
of only 27 Mpc, in the outskirts of the galaxy NGC 2272. 

Since accretion is observed most unambiguously in the 
X-rays, we used the Chandra X-ray Observatory. We ob¬ 
served SN 2012hn with Chandra for a continuous 30 ksec 
on 17 August 2013 (OBSJD=15668), 533 days after the 
explosion, which occurred on 12 March 2012. Data were 
taken in timed exposure mode and telemetered in very faint 
mode. Data redu ction and analysis we re completed using 
CI AO version 4.6 dFrusc ione et al1l2006ll and Sherpa version 
2 dFreeman et al.|[200ll l. We reprocessed the level=l event 



Figure 1. Joint upper limits for the mass of the WD and BH 
involved in the tidal detonation of SN 2012hn. Our X-ray obser¬ 
vation upper limit excludes parameter space to the upper right of 
the curves. A description of y and how these limits are calcu¬ 
lated is provided in Section |4] The horizontal jump that appears 
on the plot for rjeff = 1 is the hard/soft-state transition (the 
state transition is not crossed on the other plotted curves). The 
excluded parameter space, where the WD is swallowed whole or 
the WD produces too many heavy elements, is discussed in sec¬ 
tion (note: the boundaries are approximate). 


file using the chandra_repro script to apply the most re¬ 
cent calibration updates available (CALDB version 4.6.1.1). 

The source was not detected in the observation. To cal¬ 
culate upper limits to the source flux, we extracted source 
and background spectra using SPECEXTRACT at the loca¬ 
tion of the transient in a 6” radius circle and a 25” ra¬ 
dius annulus excluding the source region, respectively. We 
jointly fit the unbinned source and background spectra at 
0.5 - 8.0 ke V using the C-statistic, which is similar to the 
Icashi dlrrali statistic but with an approximate goodness- 
of-fit measure. Since we constrain Lx 0.02LEdd for any 
BH mass we consider , the accret i on di sk is in the canonical 
low/hard state (e.g., Maccaron3 l2003l l . in which the spec¬ 
trum is usually well-described by a power-law model. For 
the source spectral model (xsphabs x xspowerlaw), we 
set the photon index of the power-law, F = 1.6. This value 
is ty pical for accretio n in average X-ray binary populations 
fe.g. ISdl et al.l[201lll . The colu mn density, Nh was froze n 
at the galactic foreground value llDickev fc LockmanlllQOlJ l. 
The background spectrum was fit with a T = 1. 4 power law, 
consi stent with the hard X-ray background (e.g. iTozzi et al.l 
I 2 OO 6 II . We find a 3ct upper limit to the unabsorbed flux, 
Fo.s-skev < 5.1 X 10“'^® erg cm^ s“^. At a distance of 
27 Mpc, this corresponds to Lo.sskeV < 4.4 x 10®® erg s“®. 

To put this upper limit in conte xt, we adopt the an¬ 
alytical relations from ET et al.l (l2002ll for the similar case 
of the tidal disruption of a star by a supermassive BH to 
calculate joint upper limits on the mass of the WD and 
BH assumed to be involved in the disruption (given the X- 
ray luminosity upper limit above at the time of the obser- 
vation). While m ore complex numerical models exist (e.g. 
iPiran et al.ll2015l l. we opt to use the simpler analytical for¬ 
mulation to derive parameter limits, especially since the 
source is not X-ray-detected. 
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_ First, we combine equations (10) and (11) from iLi et alJ 

(l2002ti to calculate the fallback accretion rate, M, as a func¬ 
tion of time, t, fallback fraction, /, and mass of the BH, 
Mbh- We then ca lculate the acc retion efficiency, e, using 
equation (12) from iLi et al.l (l2002ll . In addition, we make a 
few more assumptions that simplify or improve our model: 

• We use the mass-radius relation for WDs, assuming 
the mean molecula r weight per free electron, /Te = 2 (e.g. 
iHansen et al. I l2004l . their equation 7.103). 

• We assume that / = 0.35, as suggested bv lROOl . 

• We take into account the accretion disk state transition 
of accreting BH binaries so that L oc rh f or L/LEdd ^ 0-02 
and L oc rh^ for L/LEdd < 0.02 (e.g. lEsin et ^ 1 1 19971: 
[m accaron i[2^). 

• To account for the fraction of light falling in the Chan¬ 
dra bandpass, we assume a bolometric correction, <5, of 0.85 
in the high, soft state and 0.1 in the low, hard state (e.g. 
iPortegies Zwart et al.l[2004fl . 


To take into account that the various scale factors are 
jointly uncertain by approximately a factor of a few, we 
define a relative effective efficiency: 


»7e// = 



I ^soft,hard \ 

V 0.85,0.1 J ^ 


( 2 ) 


We plot the joint upper limits in Mbh and Mwd for 
a range in a factor of 3 for r]eff in Fig. [T] The effects for 
different BH spins, choice of X-ray spectral model and other 
uncertainties discussed in section [U are also subsumed into 
the uncertainty in ?7e//. The fact that 0.2 < Mwd ^0.6 and 
Mbh ;S 2 X 10® Mq, which excludes large sections of param¬ 
eter space given the tidal detonation model, are discussed in 
section o and the beginning of section (2] 

From this analysis, we have placed weak constraints on 
the BH mass. It is clear that these observations are not deep 
enough or taken close enough to the time of the transient to 
place strongly constraining upper limits on the host system 
under this model. 


5 SUMMARY AND TESTS OF THE MODEL 

In this work, we have proposed that tidal detonations of low- 
mass WDs are the source of at least some of the calcium-rich 
gap transients. We have shown that the predictions of this 
model are at least qualitatively consistent with observations 
of the systems discovered thus far. Throughout this work, 
we have assessed this model using rough comparisons to test 
for plausibility; however, because of the large range of ways 
in which the tidal detonation can arise and the large uncer¬ 
tainties in key parameters in these systems, detailed quan- 
tiative comparisons are not possible at the present time and 
currently may not be computationally feasible. 

A few clear, testable predictions can be made if this 
scenario is correct: 

(i) If these transients are random flyby encounters of low- 
mass WDs with IMBHs, most of the calcium-rich gap tran¬ 
sients should be associated with low-luminosity, preferably 
dense stellar groups. However, if three-body interactions 
drive a WD already orbiting a BH into it, then these sys¬ 
tems coi^beanywhere and could include ejected systems 
(e.g. iFolevllioisI') . 


(ii) If very faint dwarf spheriodals or ultra-faint dwarfs 
are the hosts of the calcium-rich gap transients and are dis¬ 
tributed like the brighter dwarf ellipticals, then these explo¬ 
sions might occur preferentially in clusters or dense gro ups 
of galaxies (IVader fc Sandag3ll99ll : [Sabatini e t al.l l200^ . Of 
the three most recent transients presented bv lKlT two were 
located in clusters of galaxies, as was one of the two archival 
events. 

(hi) X-ray emission should be seen from the gap tran¬ 
sients, and should follow a decay law. 


Using this latter prediction, we sought to catch one of 
these calcium-rich gap transients, SN 2012hn, in the decay 
phase of the light curve. Unfortunately, we did not detect 
the source in our follow-up X-ray observations. However, 
assuming our hypothesis is correct that this transient was 
a tidal detonation event of a WD, Mbh iS 10* ^0 S’Hd 
Mwd S 0.6 Mq in this scenario. The upper limits are not 
very constraining because the effective efficiency that in¬ 
cludes the accretion efficiency, fallback efficiency, and the 
bolometric correction is highly uncertain. If this scenario is 
correct, a future X-ray observation of one of these sources 
as close as possible to the time of the transient might detect 
it and provide much stronger constraints on the masses and 
the nature of the event. 

Finally, concerning the first two points above, some ini¬ 
tial work has been completed to searc h for the hosts of a 
few of these transients without success (iLvman et al.llioi^ l. 
These recent observations rule out for a few sources, includ¬ 
ing SN 2012hn, that the transient was caused by the tidal 
detonation induced by a random flyby of an IMBH because 
the upper limits to the masses of host stellar associations are 
too low. However, the current sample o f calcium-rich gap 
transients might be heteorgeneous, and iLvman et al.l only 
placed stringent upper limits for a few sources. In addition, 
the stringent upper limits rule out most but not all globu¬ 
lar clusters and dwarf galaxies. Finally, there are numerous 
ways three-body interactions could perturb the orbits of a 
stellar system to force a WD too close to a stellar-mass BH 
or neutron star. Future detailed, multiwavelength observa¬ 
tions are clearly needed to constrain this and other various 
models for these transients. 


6 ACKNOWLEDGMENTS 

We are grateful to Giuseppe Lodato, Cole Miller, Hagai 
Perets, and Denija Crnojevic for useful discussions. This 
work has been funded through Chandra grant GO4-15061X. 


REFERENCES 

Anderson J., van der Marel R. P., 2010, ApJ, 710, 1032 
Anosova Z. P., 1990, Celestial Mechanics and Dynamical 
Astronomy, 48, 357 

Antognini J. M., Shappee B. J., Thompson T. A., Amaro- 
Seoane R, 2014, MNRAS, 439, 1079 
Antonini F., Murray N., Mikkola S., 2014, ApJ, 781, 45 
Arnett W. D., 1982, ApJ, 253, 785 

Bash F. N., Gebhardt K., Goss W. M., Vanden Bout P. A., 
2008, AJ, 135, 182 


© RAS, MNRAS 000, [THIO] 






























Calcium-rich gap transients 9 


Baumgardt H., Hut P., Makino J., McMillan S., Portegies 
Zwart S., 2003, ApJL, 582, L21 
Baumgardt H., Makino J., Ebisuzaki T., 2004, ApJ, 613, 
1143 

Bekki K., Couch W. J., Drinkwater M. J., Gregg M. D., 
2001, ApJL, 557, L39 

Belczynski K., Bulik T., Fryer C. L., Ruiter A., Valsecchi 
F., Vink J. S., Hurley J. R., 2010, ApJ, 714, 1217 
Belokurov V. et al., 2010, ApJL, 712, L103 
Bergeron P., Saffer R. A., Liebert J., 1992, ApJ, 394, 228 
Bildsten L., Shen K. J., Weinberg N. N., Nelemans G., 
2007, ApJL, 662, L95 

Binney J., Tremaine S., 2008, Galactic Dynamics: Second 
Edition. Princeton University Press 
Blecha L., Ivanova N., Kalogera V., Belczynski K., Fregeau 
J., Rasio F., 2006, ApJ, 642, 427 
Bragaglia A., Renzini A., Bergeron P., 1995, ApJ, 443, 735 
Cash W., 1979, ApJ, 228, 939 

Chatterjee S., Rasio F. A., Sills A., Glebbeek E., 2013, ApJ, 
777, 106 

Chen X., Madau P., Sesana A., Liu F. K., 2009, ApJL, 697, 
L149 

Chen X., Sesana A., Madau P., Liu F. K., 2011, ApJ, 729, 
13 

Clausen D., Eracleous M., 2011, ApJ, 726, 34 
Cseh D., Kaaret P., Corbel S., Kording E., Coriat M., 
Tzioumis A., Lanzoni B., 2010, MNRAS, 406, 1049 
Dado S., Dar A., 2013, arXiv: 1301.3333 
de Rijcke S., Buyle P., Dejonghe H., 2006, MNRAS, 368, 
L43 

Dickey J. M., Lockman F. J., 1990, ARAA, 28, 215 
Eggleton P., Kiseleva L., 1995, ApJ, 455, 640 
Esin A. A., McClintock J. E., Narayan R., 1997, ApJ, 489, 
865 

Faber S. M., Jackson R. E., 1976, ApJ, 204, 668 
Farrell S. A. et ah, 2012, ApJL, 747, L13 
Farrell S. A., Webb N. A., Barret D., Godet O., Rodrigues 
J. M., 2009, Nature, 460, 73 
Filippenko A. V., 1997, ARAA, 35, 309 
Foley R. J., 2015, arXiv: 1501.07607 
Foley R. J., Sanders N. E., Kirshner R. P., 2011, ApJ, 742, 
89 

Freeman P., Doe S., Siemiginowska A., 2001, in Society 
of Photo-Optical Instrumentation Engineers (SPIE) Con¬ 
ference Series, Vol. 4477, Society of Photo-Optical In¬ 
strumentation Engineers (SPIE) Conference Series, J.- 
L. Starck & F. D. Murtagh, ed., pp. 76-87 
Fregeau J. M., Joshi K. J., Portegies Zwart S. F., Rasio 
F. A., 2002, ApJ, 570, 171 

Fruscione A. et ah, 2006, in Society of Photo-Optical In¬ 
strumentation Engineers (SPIE) Conference Series, Vol. 
6270, Society of Photo-Optical Instrumentation Engineers 
(SPIE) Conference Series 

Fuller J., Derekas A., Borkovits T., Huber D., Bedding 
T. R., Kiss L. L., 2013, MNRAS, 429, 2425 
Gallagher, HI J. S., Wyse R. F. G., 1994, PASP, 106, 1225 
Galleti S., Bellazzini M., Federici L., Buzzoni A., Fusi Pecci 
F., 2007, A&A, 471, 127 

Gebhardt K., Rich R. M., Ho L. C., 2002, ApJL, 578, L41 
Gerssen J., van der Marel R. P., Gebhardt K., 
Guhathakurta P., Peterson R. G., Pryor C., 2002, AJ, 
124, 3270 


Gill M., Trent! M., Miller M. C., van der Marel R., Hamil¬ 
ton D., Stiavelli M., 2008, ApJ, 686, 303 
Godet O., Barret D., Webb N. A., Farrell S. A., Gehrels 
N., 2009, ApJL, 705, L109 

Grindlay J. E., Heinke C. O., Edmonds P. D., Murray S. S., 
Gool A. M., 2001, ApJL, 563, L53 
Hamers A. S., Pols O. R., Glaeys J. S. W., Nelemans G., 
2013, MNRAS, 430, 2262 

Hansen B. M. S., Kalogera V., Rasio F. A., 2003, ApJ, 586, 
1364 

Hansen C. J., Kawaler S. D., Trimble V., 2004, Stellar 
interiors : physical principles, structure, and evolution. 
Springer 

Hayasaki K., Stone N. C., Loeb A., 2015, arXiv: 1501.05207 
Hayden B. T. et ah, 2010, ApJ, 712, 350 
Heggie D. G., Rasio F. A., 1996, MNRAS, 282, 1064 
Hermes J. J. et al., 2013, MNRAS, 436, 3573 
Hodapp K. W. et al., 2004, Astronomische Nachrichten, 
325, 636 

Holley-Bockelmann K., Giiltekin K., Shoemaker D., Yunes 
N., 2008, ApJ, 686, 829 

Illingworth G., King 1. R., 1977, ApJL, 218, L109 
Ivanov P. B., Polnarev A. G., Saha P., 2005, MNRAS, 358, 
1361 

Jiang Y.-F., Stone J. M., Davis S. W., 2014, ApJ, 796, 106 
Johnston K. V., Hernquist L., Bolte M., 1996, ApJ, 465, 
278 

Johnston K. V., Spergel D. N., Hernquist L., 1995, ApJ, 
451, 598 

Kasliwal M. M., 2011, PhD thesis, California Institute of 
Technology 

Kasliwal M. M. et al., 2012, ApJ, 755, 161 
Kepler S. O., Kleinman S. J., Nitta A., Koester D., Cas- 
tanheira B. G., Giovannini O., Gosta A. F. M., Althaus 
L., 2007, MNRAS, 375, 1315 

Khan F. M., Holley-Bockelmann K., Berczik P., Just A., 
2013, ApJ, 773, 100 

Kilic M., Allende Prieto C., Brown W. R., Koester D., 2007, 
ApJ, 660, 1451 

King A. L. et al., 2012, ApJL, 746, L20 
Kozai Y., 1962, AJ, 67, 579 

Lasota J.-P., Alexander T., Dubus G., Barret D., Farrell 
S. A., Gehrels N., Godet O., Webb N. A., 2011, ApJ, 735, 
89 

Leigh N., Geller A. M., 2012, MNRAS, 425, 2369 
Leigh N. W. C., Geller A. M., 2013, MNRAS, 432, 2474 
Li G., Naoz S., Kocsis B., Loeb A., 2014, ApJ, 785, 116 
Li L.-X., Narayan R., Menou K., 2002, ApJ, 576, 753 
Lidov M. L., 1962, PLANSS, 9, 719 

Liebert J., Bergeron P., Eisenstein D., Harris H. C., Klein¬ 
man S. J., Nitta A., Krzesinski J., 2004, ApJL, 606, L147 
Lodato G., Rossi E. M., 2011, MNRAS, 410, 359 
Lu T.-N., Kong A. K. H., 2011, ApJL, 729, L25 
Lyman J. D., James P. A., Perets H. B., Anderson J. P., 
Gal-Yam A., Mazzali P., Percival S. M., 2013, MNRAS, 
434, 527 

Lyman J. D., Levan A. J., Ghurch R. P., Davies M. B., 
Tanvir N. R., 2014, MNRAS, 444, 2157 
Maccarone T. J., 2003, A&A, 409, 697 
Maccarone T. J., 2004, MNRAS, 351, 1049 
Maccarone T. J., 2005, MNRAS, 364, 971 


© RAS. MNRAS OOO.ITIITOl 


10 Sell et al. 


Maccarone T. J., Fender R. P., Tzioumis A. K., 2005, MN- 
RAS, 356, L17 

Maccarone T. J., Servillat M., 2008, MNRAS, 389, 379 
Mardling R. A., 2013, MNRAS, 435, 2187 
Marsh T. R., Dhillon V. S., Duck S. R., 1995, MNRAS, 
275, 828 

Mayer L., Mastropietro C., Wadsley J., Stadel J., Moore 
B., 2006, MNRAS, 369, 1021 
McConnell N. J., Ma C.-P., 2013, ApJ, 764, 184 
McLaughlin D. E., Anderson J., Meylan G., Gebhardt K., 
Pryor G., Minniti D., Phinney S., 2006, ApJS, 166, 249 
Metzger B. D., 2012, MNRAS, 419, 827 
Metzger B. D., Piro A. L., Quataert E., Thompson T. A., 
2009, arXiv: 0908.1127 

Michaely E., Perets H. B., 2014, ApJ, 794, 122 
Miller J. M., Raymond J., Reynolds C. S., Fabian A. C., 
Kallman T. R., Homan J., 2008, ApJ, 680, 1359 
Miller M. C., Hamilton D. P., 2002, MNRAS, 330, 232 
Miller-Jones J. G. A. et ah, 2012, ApJL, 755, LI 
Munoz R. R., Geha M., Willman B., 2010, AJ, 140, 138 
Murray N., Ghiang J., Grossman S. A., Voit G. M., 1995, 
ApJ, 451, 498 

Musielak Z. E., Quarles B., 2014, Reports on Progress in 
Physics, 77, 065901 

Naoz S., Farr W. M., Lithwick Y., Rasio F. A., Teyssandier 
J., 2013, MNRAS, 431, 2155 

Nelemans G., Yungelson L. R., Portegies Zwart S. F., Ver- 
bunt F., 2001, A&A, 365, 491 
Newell B., Da Gosta G. S., Norris J., 1976, ApJL, 208, L55 
Novikov 1. D., Thorne K. S., 1973, in Black Holes (Les 
Astres Occlus), Dewitt C., Dewitt B. S., eds., pp. 343- 
450 

Noyola E., Gebhardt K., Bergmann M., 2008, ApJ, 676, 
1008 

Ostrov P. G., Forte J. C., Geisler D., 1998, AJ, 116, 2854 
Perets H. B. et ah, 2010, Nature, 465, 322 
Phillipps S., Drinkwater M. J., Gregg M. D., Jones J. B., 
2001, ApJ, 560, 201 

Phinney E. S., 1989, in lAU Symposium, Vol. 136, The 
Genter of the Galaxy, Morris M., ed., p. 543 
Piran T., Sadowski A., Tchekhovskoy A., 2015, arXiv: 
1501:02015 

Portegies Zwart S., Boekholt T., 2014, ApJL, 785, L3 
Portegies Zwart S. F., Dewi J., Maccarone T., 2004, MN¬ 
RAS, 355, 413 

Ran A. et ah, 2009, PASP, 121, 1334 
Redmount 1. H., Rees M. J., 1989, Comments on Astro¬ 
physics, 14, 165 

Rees M. J., 1988, Nature, 333, 523 
Rhode K. L., Zepf S. E., 2004, AJ, 127, 302 
Rosswog S., Ramirez-Ruiz E., Hix W. R., 2009, ApJ, 695, 
404 

Sabatini S., Davies J., Scaramella R., Smith R., Baes M., 
Linder S. M., Roberts S., Testa V., 2003, MNRAS, 341, 
981 

Searle L., Zinn R., 1978, ApJ, 225, 357 
Sell P. H., Pooley D., Zezas A., Heinz S., Homan J., Lewin 
W. H. G., 2011, ApJ, 735, 26 
Sesana A., Haardt F., Madau P., 2006, ApJ, 651, 392 

Sesana A., Haardt F., Madau P., 2007, ApJ, 660, 546 

Sesana A., Haardt F., Madau P., 2008, ApJ, 686, 432 


Shcherbakov R. V., Pe’er A., Reynolds C. S., Haas R., Bode 
T., Laguna P., 2013, ApJ, 769, 85 
Shen K. J., Bildsten L., 2009, ApJ, 699, 1365 
Smart! S. J. et ah, 2014, ArXiv: 1411.0299 
Stabler S. W., 2010, MNRAS, 402, 1758 
Stone N., Sari R., Loeb A., 2013, MNRAS, 435, 1809 
Strader J., Chomiuk L., Maccarone T. J., Miller-Jones 
J. G. A., Seth A. G., Heinke C. O., Sivakoff G. R., 2012, 
ApJL, 750, L27 

Strader J. et ah, 2013, ApJL, 775, L6 
Strickler R. R., Cool A. M., Anderson J., Cohn H. N., Lug¬ 
ger P. M., Serenelli A. M., 2009, ApJ, 699, 40 
Sullivan M. et ah, 2011, ApJ, 732, 118 
Taylor J. M., Grindlay J. E., Edmonds P. D., Cool A. M., 
2001, ApJL, 553, L169 

Tombesi F., Melendez M., Veilleux S., Reeves J. N., 
Gonzalez-Alfonso E., Reynolds C. S., 2015, arXiv: 
1501.07664 

Tozzi P. et ah, 2006, A&A, 451, 457 
Ulvestad J. S., Greene J. E., Ho L. C., 2007, ApJL, 661, 
L151 

Vader J. P., Sandage A., 1991, ApJL, 379, LI 
Valenti S. et al., 2014, MNRAS, 437, 1519 
Valtonen M., Karttunen H., 2006, The Three-Body Prob¬ 
lem. Gambridge University Press 
Van Den Bergh S., 2003, ApJ, 590, 797 
Vennes S. et al., 2011, ApJL, 737, L16 
Verbunt F., 1987, ApJL, 312, L23 

Waldman R., Sauer D., Livne E., Perets H., Glasner A., 
Mazzali P., Truran J. W., Gal-Yam A., 2011, ApJ, 738, 
21 

Wang J., Merritt D., 2004, ApJ, 600, 149 
Wiersema K., Farrell S. A., Webb N. A., Servillat M., Mac¬ 
carone T. J., Barret D., Godet O., 2010, ApJL, 721, L102 
Willman B. et al., 2005, AJ, 129, 2692 
Woosley S. E., Kasen D., 2011, ApJ, 734, 38 
Woosley S. E., Taam R. E., Weaver T. A., 1986, ApJ, 301, 
601 

Woosley S. E., Weaver T. A., 1994, ApJ, 423, 371 
Yuan F., Kobayashi C., Schmidt B. P., Podsiadlowski P., 
Sim S. A., Scalzo R. A., 2013, MNRAS, 432, 1680 
Zaritsky D., Gonzalez A. H., Zabludoff A. L, 2006, ApJ, 
638, 725 

Zinn R., 1985, ApJ, 293, 424 

This paper has been typeset from a T^X/ DTf^X file prepared 
by the author. 


© RAS. MNRAS OPO.IHflOl 


